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Abstract 

We deal with the 3D inviscid Leray-a model. The well posedness 
for this problem is not known; by adding a random perturbation we 
prove that there exists a unique (in law) global solution. The random 
forcing term formally preserves conservation of energy. The result 
holds for initial velocity of finite energy and the solution has finite 
energy a.s.. These results are easily extended to the 2D case. 
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1 Introduction 

The motion of incompressible fluids is described by the Navier-Stokes equa- 
tions 

— -uAv + {v V)v + yp = f 
div V = 

for viscous fluids, or by the Euler equations 

div V = 



*Universita di Padova, Dipartimento di Matematica, Padova, Italy, bar- 
bato@niath.unipd.it 

^University of Wyoming, Department of Mathematics, Dept. 3036, 1000 East Univer- 
sity Avenue, Laramie WY 82071, United States, bessaih@uwyo.edu 

^Universita di Pavia, Dipartimento di Matematica, via Ferrata 1, 27100 Pavia, Italy, 
benedetta.ferrario@unipv.it 



1 



for inviscid fluids. 

The unknown are the velocity field v = v{t,x) and the pressure field 
p = p{t,x); / is a given external force and u > is the viscosity that 
corresponds to the inverse of the Reynolds number Re. When the fluid 
moves in a bounded domain, suitable boundary conditions are associated to 
these equations, respectively, the no-slip and slip conditions. 

The above two systems have a quite different behavior; for instance, when 
/ = system ([1]) is dissipative while system ([2]) is conservative. It is well 
known, since the seminal work of Leray, that for initial velocity of finite 
energy the 3D Navier-Stokes system ([1]) has a global weak solution but its 
uniqueness is still an open problem. However, for the 3D Euler system ([2]) 
neither the global existence nor the uniqueness of global solutions are known 
when the initial velocity is of finite energy (we refer to the review paper [4] 
on this topic). 

We recall that to prove the existence of solutions to ([T]) in M*^, d = 2,3, 
Leray [23] considered the following regularization for a > 

dv°' 

— - z/At;° + {u" ■ V)w° + yp = f 
div v°' = 

where Ga is a smoothing kernel such that u°' — > v^, in some sense, as a — 0. 
In particular, system (|3]) converges to the Navier-Stokes system ([1]) as a — 0. 

In [8J, a special smoothing kernel was considered, namely, the Green 
function associated with the operator 1 — aA, 

u" = Ga*v'' = {l-aA)-^v'' (4) 

for a > 0. This kernel works as a kind of filter with width a and the param- 
eter a reflects a sub-grid length scale in the model. This model was inspired 
by the Navier-Stokes a model (also known as the Camassa-Holm system or 
Lagrangian averaged Navier-Stokes a equations) of turbulence, see [6l Ul [E] 
and the references therein. Moreover, it has been demonstrated analytically 
and computationally that the Navier-Stokes a model is a powerful tool in 
the study of turbulence, see [12] and the reference therein. Along the same 
lines, it is worth mentioning that other a models, such as the Clark-a model 
[5] and the Navier-Stokes Voigt equations [20], have been used as a sub-grid 
scale models of turbulence. 

In this paper, we are interested in a stochastic version of system with 
1/ = and regularization given by (jl]), that is the following stochastic Leray-a 
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model of Euler equations 

dv°' + [(m" ■ S/)v°' + Vp] dt = ((cr o dW) ■ V)w" 
y" = (l-aA)u° (5) 
div f ° = 

Here is a Brownian motion (in time) and odW referes to the Stratonovitch 
differential; the parameter a is positive. When a = the first equation of 
([S]) reduces to the stochastic Euler equation. 

The well posedness of weak solutions of the deterministic system (jSD (c" = 
0) is not known. In particular, when the initial velocity has finite energy, 
existence of global weak solutions can be proven if a > 0, see the Appendix. 
However, the uniqueness is not known for both d = 2 and d = 3. If a = 
and (7 = 0, global existence of solutions are known for initial velocity of finite 
energy and enstrophy for d = 2 while it is open for d = 3; their uniqueness 
is an open problem for = 2 and d = 3 (see, e.g., [Ij). 

Let us point out that the analysis of the deterministic Leray-a Euler 
equations in 3D, that is system (5) with a = 0, is more difficult than for 
the other approximations models for 3D inviscid fluids (i.e., Camassa-Holm, 
Clark and Voigt models). Indeed, for these other models there is formal 
conservation of the sum J \v°'\'^dx + a J \'Vv°'\'^dx {a > 0), whereas the model 
we are interested in has only formal conservation of / |f From this point 

of view the deterministic Leray-a model for 3D Euler equations considered in 
this paper is closer to the 3D Euler equations than the Voigt, Camassa-Holm 
and Clark models which regularize much more the original Euler equations. 

When adding an appropriate stochastic perturbation, we will prove that 
system ([5]) has a unique global solution (in law) when the initial condition is of 
finite energy. We will prove existence and uniqueness (in law) of solutions by 
means of Girsanov formula. The multiplicative noise in ([5]) formally preserves 
the conservation of energy (see Section 3 for the details). It is crucial to 
choose the random perturbation in to be written in the Stratonovitch 
form in such a way that formally the energy of the vector field is conserved. 

All our results are stated for a three dimensional spatial domain (a box 
[0, 27r]^, assuming periodic boundary conditions), but our proofs can be easily 
adapted to the two dimensional case. It would be interesting to study the 
behavior of the process when a and/or a converge to zero. This is the 
subject of future research. 

In the past few years, there has been a huge effort to tackle the problem 
of using a similar noise in order to improve the qualitative properties of some 
non linear equations. In particular uniqueness of the stochastic equation 
has been proved either when uniqueness is not known in the deterministic 
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setting or with weaker assumptions than in the deterministic setting; for 
these results, see PQ 121 El [HI [131 [IH [15] and the references therein. We refer 
to [HI [16] for the 2D Euler equations, and to [101 [E] for some analysis on the 
3D Navier-Stokes equations. For an overview of the problems and methods, 
we refer to [12]. 

As far as the content of this paper, in Section 2 we will introduce our 
functional setting and spaces. Section 3 will be devoted to the description 
of the stochastic Leray-a model in Fourier components. We will write the 
model in both the Stratonovitch and Ito forms. Section 4 will focus on the 
linear model: global existence and uniqueness of strong (in the probabilistic 
sense) solutions will be proved. The uniqueness proof is based on the study of 
a new linear problem constructed by means of the covariance matrix A^. The 
nonlinear model will be studied in Section 5, where we prove the existence and 
uniqueness of solutions (in law) by means of a Girsanov formula. In Section 6, 
our results will be stated for the stochastic partial differential equation ([5|). 
To make our paper self-contained, in the Appendix we will give the proof 
of global existence of weak solutions for the deterministic Leray-a model of 
Euler equations. 

2 Functional setting 

Let the spatial domain be a torus T, i.e. x G and periodic boundary 
conditions on the cube [0, 27i]^ are assumed. Notice that if v is a solution, 
then also f + c (c G M) is a solution. Therefore, we consider mean zero 
velocity vectors, i.e. /[q 27r]3 "^(^5 ^-^ ~ ^■ 

We fix notations. Given a complex number 2; G C, we denote, respec- 
tively, by and its real and imaginary part; hence, z = ^z — i'^z and 
the product of two complex numbers z and w is wz = (JHw^z — '^w'^z) + 
i{^w^z + '^w^z). 

Morever, let x G be represented (^)). For x,y e C^, 

we set {x,y) = X]j=i ^^"^ II^P — {x,^)- This defines, as a particular 

case, also the scalar product and the norm in M.^. 

For each k G Z^, let ek{x) := e^^^'''\ x G M^. The family {ek}k£Z^ is a 
complete orthogonal basis for the space L^{T, C). 

In this Section, we write the deterministic Leray-a model of the Euler 
system ([5|) 
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in Fourier components; this will be given in ([H]). In the next Section, the 
stochastic forcing term will be introduced. 

Assume v(t, ■) and u(t, ■) are in {L'^(T, C))^; then 

v(t,x) = ^ Vkit)ek{x) and u{t,x) = ^ Uk{t)ek{x) 

We have = and Uq = 0, since v has mean value zero. Moreover, since v 
and u are real valued and e_fc(a;) = ek{x), we have 



We set \\v{t,.)\\l = j:k\\Mt)r- 
We substitute in ([6j)2 and get 

^Vkit)ek{x) = (1 - aA) ^Mfc(t)efc(x). 

k k 

From now on, we will drop the index a in the unknowns for simplicity. 
Since Aefc(x) = — ||A;|pefc(a;) then ffc(t) = (1 + a;||A;|p) Uk(t) and 

'^Vkit)ek{x) = ^ (l + a\\kf) Uk{t)ek{x). 

k k 

From equation we get the incompressibility condition 

{vk{t),k) = \/k eZ^, t eR. 
Finally, using ([6])i we obtain 

^^Mt)ekix) 

k 

^Uh{t)eh{x) ■ V I ^ek'{x)vk'{t) - \/p{t,x) 

\ h / k' 

-^e,(x) (u^^\t)^ + u^^\t)g^ + u^^\t)g^') ekix)vkit) - yp{t,x) 

h,k' 

-J^^hix) (4')(t)^A;'« +nf (t)zA;'(2) +43)(^),^/(3)j ek'{x)vk'{t) - Vp{t,x) 

h,k' 

-i^eh+k'{x){uh{t),k')vk'{t) - Vp(t,x), 



h,k' 
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that is 



di 



where 



Pk(v) :=v 



{k,k) 



k 



(7) 



is the projection onto the space orthogonal to k. 

Summing up, since {vh, k) = {vh, k — h), we obtain the system (P) written in 
Fourier components 

dvk. . {vh{t),k) 



{vkit),k) = 
for any k. 

We notice that, given a > 0, if Ylik < then the series in the 

r.h.s. of ([H])i is convergent and for any k we have 



<\\vmm\- 

System ([8]) enjoys an important property: formally the energy E{t) 



\ Ylik is conserved under the dynamics given by ([8]). Indeed, 



d „ 



{vh{t),k) 



+ a\\h\ 



[Pk{Vk-h{t)),Vk{t)) 



Summing over all components, we formally obtain 



dE 
~dt 



k h ^ 



+ a 



{Vk-h{t),Vk{t)) 



which vanishes, since the sum contains terms which cancel each other ac- 
cording to the following equality 



———{Vk'-h',Vk'} = ——-r7r^{Vk-h,Vk} 

1 + a /i r 1 + a /i r 



(9) 



for h' = —h and k' = k — h. 

Let us finally notice that conservation of energy formally holds also for 
a = 0. 
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3 The stochastic Leray-a model in Fourier 
components 

We are interested in a stochastic equation obtained from ([8]) by adding a 
random forcing term in such a way that energy is formally conserved. To 
this end we consider the system of Stratonovich equations 

h&Z^ " " 

(10) 

where {Wh}h&^ is a family of independent C'^-valued Brownian motions on a 
filtered probability space {fl, {J-'t)t>o, P), except for {Whit), h) = 0, W-h{t) = 
Wh{t). According to the properties of Stratonovich integral (see [19]) we have 
formally that dE{t) = 0. The computations are similar to the previous ones, 
using and 



{Yk-h,Yk){crodW,h,k) = {Yk,Yk.H){aodWh,k). 

Let us make precise the Stratonovich formulation of system (ITU]) in order 
to write it in terms of Ito integrals. Indeed, the Stratonovitch formulation 
gives insights on the behaviour of the system, but computations will be done 
on the Ito formulation. 

Set 

J ■= {{ki, k2, ks) G : fci > or (fci = 0, fcg > 0) or (fci = 0, /cs = 0, k^ > 0)}. 

Let {Wlj}h£j be a family of independent C^-valued standard Brownian mo- 
tions. Define for any h & J 

Wh:=Ph{Wi), W.h.= Wh (11) 

Therefore (Wh, h) = 0, W^h = Wh for any h E , h 0. Hence it is enough 
to give the family {Wl^}h£j in order to define the stochastic part of system 

Next, define 

Wh,k{t) ■■= {Wh{t),k) (12) 
then each Wh,k is a C-valued Brownian motion, whose real and imaginary 

part are independent. Since WhA^) = {Whit),k) = {Ph{Wi{t)),k) = {W[Xt) , Ph{k)) , 
then 

Var[mVhM] = Var[^WhM] = \\Phik)f = WkW" sin" (9) 
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where 6 is the angle between h and k. Now, setting 

WhA\Ph{k)\\ ■.= Wn,k (13) 

we have defined standard real Brownian motions W//h,k and ^Wh,k- 
Set 

Therefore (HOl) is 



(n(t),A:) 

] 

(yh\\Ph{k)\\Pk{Yk-h{t)) o dWh^k{t) 



(15) 



I 



Indeed, for the Stratonovich integral we have 

t 

Pk{Yk-h{s)){ahodWh{s),k) 

= ahPk [ Yk^h{s) o dWh,k{s) 
Jo 

= ah\\Ph{k)\\Pk I Yk-h{s) o dWh,k{s) 
Jo 

= a4Ph{k)\\Pk [ ^Yk-his) o dmVh,k{s) - (rh\\Ph{k)\\Pk I '^Yk-nis) o d'^WhAs) 



+ ta,,\\Pf,{k)\\P, / m-h{s)od^Wh,k{s) + ta4Ph{k)\\Pk / ^Yk-his)od^W,,^k{s). 
Jo Jo 

(16) 

We have the corresponding Ito formulation. 

Theorem 1. Let {i^fcjfcgzs fc^g ^ sequence of continuous and adapted pro- 
cesses defined on a given filtered probability space such that J2k ll^fc(^) P < oo 
a.s.. If the sequence solves the following system 

dY,{t) = - ^ V TT^%?|^fc(n-/.(t))rft 

- I J2 ^h\\Ph{k)\\PkiYk^h{t))dWh4t) (17) 

- J2 4\\PhikWPk{Pk-h{Ykmdt, \/k 

then it solves system f|T5l) . 
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Proof. We are going to prove that when the Ito integral in the r.h.s. is written 
as a Stratonovich integral, we get f[T^ . The corrective term appearing in this 
transformation comes from the quadratic variation [(Th\\Ph{k)\\Pk(Yk-h) , Wh,k] 
(see [igj). 

Let us work on the real and imaginary part; this makes the proof long 
but clear. Let k' = k — h; from (HM we have 

ah\\Ph{k)\\Pki^Ykit)) = ah\\Ph{k)\\Pk{^Yk>{0)) + [ ds 

Jo 

+ [ ^h\\Ph{k)\\Pk{ah'\\Phik')\\Pki'^Yk>_His)))odm^^^ 

fT;.||P^(A;)||Pfc(53n,(t)) = a^||P^(A;)||Pfe(53n.(0))+ / ds 

Jo 

-Y. I (yh\\Phm\Pk{oh.\\Ph\k')\\Pk\^Yk:-h\s)))od^w^^ 

+ Y. I ^h,||P/^(A;)||PfeKHI^/.K^0ll^feK^>n'-/^K^)))°^^^>^^/^',fc'(5)- 

Bearing in mind that 

Wh',k-h is independent of Wh,k if h' ^ —h and h' ^ h 

and ^ ^ 

Wh',k-h = Wh,k if h' = -h. 

Wh',k-h = Wh,k if h' = h 

we are reduced to take into account only the terms with h' = —h and h' = h. 
Moreover, we use that Ph{k) = P~h{k ~ h) = Ph{k — h). Therefore the 
corrective term for cr/i||P/i(fc)||Pfc(3?Ffc_/i(s)) o d'^Wh,k{s) is 

-l^l\\Ph{k)fPk {Pk-h {mis))) + ^al\\P^,{k)fPk {Pk-h im-2his))) ; 

that for ah\\Ph{k)\\Pk{'iYk-H{s)) o dmVh,k{s) is 

-^iWPHm'Pk [Pk-h (mis))) - ^ai\\p,{k)rPk [Pk^H mk-2H{s))) ■ 
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that for -a,,\\PhmPk{^yk-h{s)) o dmV^,^k{s) is 

-^al\\P,ik)fPk {Pk-h C^Y.is))) - ^al\\P,{k)fPk {Pk-h {^Y,^2h{s))) ; 
that for ah\\Phik)\\PkiQYk^h{s)) o dQWh,kis) is 

-^al\\P,ik)fP, (P,_, + ^al\\P,ik)fP, (Pfc_, i'^Yk-2his))) . 

Summing up all the contributions, we get the expression given in the Propo- 
sition. □ 
The aim of this paper is to study the stochastic system (fTTI) with initial 
data of finite energy. 

4 The linear model 

Let us consider the linear system obtained by neglecting the nonlinear terms 
in (dZD: 

' dYkit) = -iEh&^^h\\Phik)\\PkiYk-hit))dBf,4t) 
-EHe^^^l\\Ph{k)rPk{Pk-h{Y,mdt 

, {Ykit),k) = ^^g^ 

Y^k{t) = W) 

. Yk{0) = Vk 

for each k ^ ~^ . Here, {B^ k} is a family of C-valued Brownian motions 
obtained from a family of independent C^- valued standard Brownian motions 
{-B^j/igj defined on a filtered probability space (fi, {J^t\t>Oi Q) with the same 
procedure presented in f|TT|) - f|T3|) . 

In the next section, we will see how Girsanov transform allows to pass 
from the linear to the nonlinear system. 

Notice that if E'^ \\Y{t)\\i2dt < oo, then the terms in the r.h.s. of f lTS]) i 
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are well defined. Indeed, for the Ito integrals we use that 
E'^WY, (^h\\Ph{k)\\ f PuiYu-h{s)dBhAs)f 



alA;fE« f \\Y{s)\\lds 
Jo 



and for the deterministic integrals (Q-a.s.) 

< r II E ^mm'PkiPk-.iY.ismds 







with E/.ez3 (^l = Y.h&? (i+4|fe|P)^ < +'^- 

We are interested in the stochastic system f|T8|) with deterministic initial 
data y = {yk}k of finite energy. We shall deal with strong solution in the 
probabilistic sense, that is the filtered probability space {Q, {J^t}t>o,Q) and 
the Brownian motions {-B^}/igj are given a priori. We shall prove existence 
and uniqueness of solutions of the following type. 

Definition 1. Given y Et^ , an energy controlled strong solution for system 
f|T8|) is a family of continuous and adapted C^-valued stochastic processes 
{Yk}kez'i such that for all t > 

' Ffc(t) = yk - iY.h&^ ^h\\Phik)\\ Jl^ Pk{Yk-his))dBh.kis) 
- E.ez3 ^l\\Pk{k)r /o Pk{Pk^HiYk{s)))ds 

{Yk{t),k) = 
I Y_k{t) = W) 
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Q-a.s. for all k, and for any t > 



Y.\\Yu{t)r<Y.\\yur Q-a.s. 



4.1 Existence of a strong solution 

Theorem 2. For any initial data of finite energy, there exists an energy 
controlled strong solution to system f|T8l) . 

Proof. We consider the finite dimensional system associated to flTSl) : for each 
integer > this is obtained by neglecting the components of index higher 
than in such a way that the energy is conserved. Set = {/i G : < 
||/i||<A^, 0<||/i — A;||< A^}. Then the linear Galerkin system is 

dYkit) = -tEh^r%^h\\Ph{k)\\Pk{Yk-hit))dBf,^kit) 



{Yk{t),k) = 



(19) 



r_,(t) = Y{t)k 

, n(0) = yk 

for each k G Z'^, < ||A;|| < A^. We consider inital data Y^{0) = y^ obtained 
from the inital data y of the full system ( 1181) by putting to the components 
y^ with > A^. 

By linearity the Galerkin system has a unique global strong solution 
Y^ = {^fc^}o<||fc||<Af- Each component is a continuous and adapted process. 
Moreover, energy is conserved, that is for any t > 



\M<N 

To prove it, again we use the properties of the family {Bh.k} of Brownian 
motions and the projector defined in ([7]). 

Therefore, for an inital data of finite energy we have for any t > 

||F^(t)||p = ||l/^||p<||l/||,2 Q-a.s. (20) 

This implies that for any p G [1, C)o) we have 

E« r WY^'m^.dt = Ib^llf, < Ibllf, VAT. (21) 

^0 



12 



Therefore the sequence {Y'^}^ is a bounded sequence in L^^Q x [0,T];/^) 
for any 1 < p < oo. This imphes that there exists a sequence {Y^^}'^-^ and 
a process Y E L°°{nx [0,T]; Z^) such that 



hm = y weakly in L^i^ x [0, T];l^) for p < cx3 

and 



hm Y^' = Y ★-weakly in L^{n x [0,T];P). 

i—^oo 

In particular 



lim Y^' = Yk weakly in L^(fi x [0,T]). 



2— >-00 



Now we consider the convergence of the integrals in the r.h.s. of (1191) i . 
The Ito integral, considered as a linear operator, is strongly continuous in 
L^(fi X [0,T]); hence it is weakly continuous (see .e.g. [22l[2l]). This implies 
that each stochastic integral converges weakly: 

lim / Pk{Y^_^{s))dBu,k{s) = [ Pu{Yu-h{s))dBh.k{s) weakly in (fix [0, T]) 



On the other side, using the independence of the Ito integrals and the Ito 
isometry we have 



\\h\\<N 



\\h\\<N 
\\h\\<N 

<a^iiA:f E E« fm\{s)rds 

„U„^AT Jo 



\\h\\<N 
,2|| Lll2||„.||2 



<^ ll^ll IIz/IIp bydSD. 
Hence 



lim E / ^h\\Ph{k)\\Pk{Yjf_j,{s))dBt,^k{s) = y^ [ a,||P,(A;)||Pfc(n_,(s))d5,,fc(s) 
weakly in L^{n x [0,T]). 
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For the deterministic integral, it is an easy computation to identify the 
hmit: 



Pk{Pk-h{Y,''{s)))ds ^ / Pk{Pk-h{Yk{s)))ds weakly in L\Qx [0, T]). 

Jo 

For the limit Y, we have for any t > 

||^WII/2 < \\y\\i2 Q - a.s. 



Therefore Y is an energy controlled strong solution. 
Moreover, using again the estimate 



/ \\Y{t)\\ldt < OO 



in a classical way we obtain that the process given by the infinite sum of 
Ito integrals '^h'^fi'\\Ph{k)\\ Pk{Yk-h{s))dBh,k{s) has a continuous modifi- 
cation. Hence, we conclude that the process Y has a continuous modification. 
□ 



4.2 The covariance matrices 

In Theorem [2] we have proved existence of energy controlled solutions Y = 
{^fc}fcez3 of ffTSl) : now we want to show their uniqueness. The idea is to study 
the time evolution of the covariance matrices {Ak}k^x3, defined as follows: 



Ai'^^'it) = E« ^Yl''\t)^Yi''\t) + QY^''\t)'^Y^''\t) 



ji,j2 = 1,2,3 



We collect the properties of A^. Since {Yf:(t),k) = for any t and k, then 
k is an eigenvector for Ak{t) corresponding to the eigenvalue. Ak{t) is a 
symmetric and semi-positive definite matrix; therefore the trace of Ak{t) is 
non negative. Moreover, we have 

J2 Tr(A,(t)) < \\y\\l (22) 

for any t > 0. Finally, PkAk{t)Pk = Ak{t), where Pk is the real matrix pre- 
viously defined in ([7]), which is symmetric semi-positive definite; P^ has the 
eigenvalue with eigenvector k and the eigenvalue 1 of double multiplicity. 

Bearing in mind fll8l) and the properties of the Brownian motions Bh^k, 
with some long but easy computations we get that each Ak fulfils a linear 
equation. 



14 



Proposition 3. For each k ^ 0, A^. fulfils the differential equation 



dAk 



(t) = - ^ al\\Ph{k)fPkPk-hAk{t) 



dt 



(23) 



+ 2Y,^l\\Ph{k)fPkAk-h{t)Pk 



This Proposition shows the non trivial fact that the covariance matrices 
satisfy a closed differential system. 

4.3 Pathwise uniqueness 

We prove pathwise uniqueness for system f|T8l) . 

Theorem 4. There exists at most one energy controlled strong solution to 
system (1181) . that is given two energy controlled strong solutions Y[i] and Y[2] 
to system ffTSjl defined on the same probability space (fi, {J^t\t>o, Q) and with 
the same initial data y E l"^ and Brownian motions, we have for any t > 



The idea of the proof is to take the difference Y = Y^i] — Y[2] ; by linearity Y 
solves (ITSj) but with initial data 1^(0) = 0. Let {Ak}k&3 be the covariance 
matrices of Y; these matrices satisfy (1231) with zero initial condition and 
regularity (^^. Thus, the uniqueness problem for (llSp is transformed in the 
easier uniqueness problem for the deterministic system fl23l) . Indeed, in order 
to show that for any t > we have Y{t) = Q-a.s. it is enough to prove 
that system fl23|) with zero initial condition has the unique solution which 
vanishes, i.e. for any k ^ 0, given ^^(0) = we have Ak{t) = for all t > 0. 
For any T > define 



Yii]{t)=Y[,]{t) 



Q — a.s. 



Proof. Define 




Each tensor Bk enjoys the same properties of A^. Since 




(24) 



15 



then 

J]Tr(i?,)<4||y||2T; 

thus 

hm Tr(5fc) = 0. (25) 
Writing equation f l23|) in the integral form we have 
AuiT) = ^l\\Ph{kW [2PkBk-hPk - PkPk-hBk - BkPk-hPk] ■ (26) 

If Bh = for all h, then Ak{T) = and the proof is completed. 

We proceed by contradiction. Suppose that there exists ki such that Bk^ 
does not vanish; then the maximal eigenvalue A* of B^-^ is positive. Starting 
from Bk^ and we construct a sequence of Bk„ and A*'s (with A* the 
maximal eigenvalue of Bk„) such that {A* is a stricly increasing sequence. 
Therefore lim„^oo A* > 0. On the other hand, each Bk is semipositive definite 
and therefore Tr(_Bfc^) > A*. It follows that 

lim Tr(5fcJ > 

n—^oo 

which is impossible because of fl2^ . 

Therefore we are left to prove that given G N such that Bk„ has 
maximal eigenvalue A* > 0, then there exists kn+i G N such that -Bfe^^^ has 
maximal eigenvalue A*_,_^ > A* > 0. 

Let (pn be the eigenvector corresponding to A*. Therefore 

= A>„, (27) 

Moreover 0„ is orthogonal to since they are eigenvectors corresponding 
to different eigenvalues; therefore 

Pkjn = K (28) 

From fl26|l we have 

< Afc„(T)0„) = ^^||n(A;)f (2(0„,,Pfc„i?,_;,P,„0„) 

— {<Pri,Pk„Pk„-hBkn(t)n) 

- {(Pn, Bk^Pk^^hPkM) 

Using that each B^ is symmetric and (1271) - (!28l) . we get 

< 2 5^ a2||P,(A;)f [(0„,i?,_^0„) - A:(0„,Pfe„_^0„)] 
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For fixed n, we have that {(pn, Bk„^h4>n) tends to as \\h\\ — )■ oo, because of 
()25p . Therefore, some addends in the sum are negative; if the sum must be 
non negative there must exist at least one addend positive, i.e. 

: (0„, - A:(0,, Pk„.-M > 

Set kn+i = kn — h. Then 

Setting ipn = Pk„+i4>n and using that Bk„^-^ = Pk„+^Bk„+^Pk^^^ we get 

which imphes that the maximal eigenvalue A*_,_^ of -Bfc^+i is larger than A*. 
□ 



5 The nonlinear model 

Consider the nonlinear system in the Ito form 

- 1 (yh\\Ph{k)\\Pk{Yk-h{t))dBh,k{t) 

- J2 cTl\\PH{k)fPk{Pk-h{Ykmdt (29) 

{Yk{t),k)=0 
Y.k{t) = Yk{t) 

. n(o) = Vk 

Starting from the solution of the linear system f lTSjl we construct a solution 
to this nonlinear system by means of Girsanov transform. We shall deal with 
solutions on any fixed finite time interval [0,T]. 

Definition 2. Given y & P , a weak solution of equation (129]) in P is a 

filtered probability space {yL,{J^t}te[Q,T],P), sequence of independent C^- 
valued Brownian motions W = {Wl^}hej on {Q, {J^t}telo,T], P) and an t^- 
valued stochastic process Y := (Yk)kez'-^ on (fi, {J^t}teYd,T]i P), with continuous 
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and adapted components Y^. such that for all t G [0, T] 

I ci i- ~r Oi\ \ tin 



— i 



(n(t),A;) = 



[ r_fc(t) = n(t) 

P-a.s. for each k G Z'^. VKe denote this solution by {{Q, {J^t}t€io,T], P), Y, W). 

Moreover, it is called an energy controlled weak solution if for all t G [0, T] 
this solution satisfies 

As usual, the Brownian motions Wh,k are constructed by the family 
{Wj^jheJ (see (IIID-(II3D). 

First we present the Girsanov result. Let Y = {Yh}^ be the strong energy 
controlled solution of (IT5I) . Define 

L{t) = -y2 fYh{s)dBh{s) 

Then L is a martingale and its quadratic variation [L, L\ is well defined and 
given by 

[L,L]{t) = ^^f^Y.^n{sfds (30) 



h&j 



We have 



Proposition 5. LetY he the strong solution of system ffTSj) with the family of 
independent 'C^ -valued standard Brownian motions {B'^}h&j on (fi, { J-i}t>o, Q) ■ 
Then ^ 

Wi{t) = B',Xt) -- I Yh{s)ds, heJ,0<t<T, (31) 
o- Jo 
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defines a family of independent -valued standard Brownian motions on 
{Q, {J^t}o<t<T, P) with the measure P, defined on {Q, J^t), which is absolutely 
continuous with respect to the measure Q and 



holds true. Therefore, Girsanov transform now apphes in a classical way (see, 



We point out that ^ > 0, Q-a.s.; therefore also ^ is well defined. Thus 
the measures P and Q are equivalent (i.e. absolutely continuous to each 
other) . 

Our main result is 

Theorem 6. For any initial data of finite energy, system fl29|) has an energy 
controlled weak solution. Moreover, this solution is unique in law. 

Proof. As far as the existence is concerned, we have that (ITSIl has a unique 
strong solution Y defined on (fi, Q) and satisfying (ITSl) . Using the 
Girsanov transform fl3T|) . we get that ((fi, {J-i}, -P), Y, W) is a weak solution 
of fl29|) . Moreover, the measure P is equivalent to the measure Q; then 
||^(^)||p < II y 11/2 P- and Q-a.s. This means that this weak solution is an 
energy controlled solution. 

As far as the uniqueness is concerned, if there were two different weak 
solutions of the nonlinear system (l29l) . then each of them would give rise to 
a weak solution of the linear system (|T8|) : these are obtained starting from 
((fi, {J^t}, P), Y, W) and getting ((^], { J^J, Q), Y, B') via Girsanov transform. 
On the other side, the pathwise uniqueness for the linear system (ITSI) implies 
the weak uniqueness; this comes from Yamada-Watanabe theorem, which is 
usually known for finite dimensional systems but whose validity holds also 
in the infinite dimensional setting as soon as the Ito stochastic integrals are 
well defined. Now, using the absolute continuity of P and Q, we deduce that 
the nonlinear system fl29|) has a unique solution (in law). □ 

Remark 7. i) The proof shows that our technique can he applied for any a > 
to more general models, that is we can deal with a noise defined by means 
of cfh = (i-|-aj[fe||2)p one? with the smoothing term given by = (1 — Q;A)~^f", 
for any p > A/3. 

a) In the 2-dimensional case, all our computations can be extended for p > 




(32) 



e.g. mm)- 



□ 



1/2. 



19 



6 The formulation in SPDE 

The stochastic model considered so far in Fourier components can be written 
as a stochastic partial differential equation, as follows 

dv + {u- V)v dt +Vpdt = E,6Z3 EU ^^^^^ ° d^h^ 
V = {1 — aA)u /oo\ 
divv = ^ ' 

viO) = vo 

For simplicilty we have dropped out the index a in the unknowns. 

The first equation of system fl5^ can be written in a more compact form 

as 

dv + Y iTT-T u^'^ dt +ypdt = y o dW^^\ (34) 
i=i i=i 

Here the random field W is given as Wit^x) := J2hez^ o'heh{x)Wh(t). 
More precisely, the Wiener process W has the following form 

^ cos((/i, xmW.jt) - sin((/^, 

^(^'^) = 2-l^ Tmnp ' ^^^^ 



where {VT/ij/ieJ is a family of independent C^-valued Brownian motions on a 
filtered probability space {Q, {J-'t)t>o, P), such that (Whit), h) = 0. 

Let us denote by H and V the subspaces of (L^(T))^ and (.^/^(T))^ re- 
spectively, given by vectors fields divergence free and periodic: 

H = {v E [L'^(T)Y , V ■ f = 0, / t>(x)(ix = 0, f ■ n periodic on T} 

Jt 

V = {v e H -.v e[H\T)f, M periodic on T} 

where n is the unit normal to the boundary of the spatial domain. 

Moreover, identifying H with its dual H' we get the Gelfand triple (V, H, V') 

V C H ^ H' cV. 

The norms are inheritaed from the spaces (L^(T))^ and (iJ^(T))^. 

Definition 3. Given vq E H , a weak solution of (l33l) in H is a filtered 
probability space {Q, {J^t}, P), o sequence of independent -valued Brownian 
motions {Wh}h on (fi, {J-'t},P) and an H-valued continuous and adapted 
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stochastic process v on {Q, {J^t},P), such that 



v{t, x), (f){x))dx — / / {{u{s, x) ■ V)(f){x),v{s, x))ds dx 
T Jo Jt 

{vo{x), (j){x))dx 

3 



X^^'^Xl / ( / ^hix){-^^{x),v{s,x))dx^ odWlf\s), P-a.s. 

(36) 



for each t G [0, T] and for all test functions : — )• M^, periodic on T, diver- 
gence free and of class . We denote this solution by {Ft}, P), v, W). 
Moreover, it is called an energy controlled weak solution if for all t > this 
solution satisfies 

\\v(t, ■)\\h < \\vo\\h P - a.s. 



This weak formulation corresponds to the stochastic equation flMj) . In- 
deed, for a more regular solution v, by integration by parts in f l5B]) one gets 
fl5^ . This is a classical result for the Euler equation and the stochastic part 
uses the properties of the Brownian motions. Therefore we have the following 
result. 

Proposition 8. The equality v{t,x) = Ylh^i^ Yh{t)^h{,x) relates the solutions 
of the stochastic PDE f l5B]) and of the stochastic Fourier system ffT^ (or 

Finally, we can reformulate our result for the SPDE: 

Theorem 9. For any initial velocity of finite energy, equation (IMI) has an 
energy controlled weak solution. Moreover, this solution is unique in law. 



7 Appendix 

In this section, we present a proof of global existence of a weak solution for the 
deterministic system ([6]) with initial velocity of finite energy; here weak has 
to be understood in the sense of PDE's. No written proof has been found in 
the published literature, whereas there are results of local (in time) existence 
and uniqueness for very regular initial velocity; however, uniqueness of weak 
solutions is an open problem. Anyway, Edriss Titi has presented a proof of 
global existence of weak solutions, as a private communication. 
For simplicity, let us drop the index a. 



21 



Theorem 10. (Due to E. S. Titi) Let vq E H and T > 0. Then there exists 
a global weak solution v to the system such that 



and 



veL°^i[0,T];H)nCi[0,T];r) 



v{t, x), (f){x))dx — / / {{u{s, x) ■ 'V)(t){x),v{s, x))ds dx 



JT 



(37) 



(t'o(x), (j){x))dx 



for each t G [0, T] and for all test functions : — )■ M^, periodic with period 
box T, divergence free and of class . 

Proof. Let P/v be the finite dimensional projector that is = P^v means 
v^{x) = '^\\^\^j^{fjv{x)eh{x)dx)eh{x). Then, we get the following finite- 
dimensional system system 

^ + p„|(«".vK| = o 

+ aAu^ = 

From (ESDi, we get j-^\\v^ it)\\jj = for any and t; then \\v^it)\\jj < hoWn- 
Thus 

sup ||v^||Loo(o,r;H) < oo; (39) 

iV 

hence using (l38|l o 

,N 



supIIm ||loo(o,T;(h2(t))3) < oo- (40) 

N 



Using again equation 



dv^jt) 

dt iuii„=i 



v,= sup |((«^(t)-V)t;^(t),P^</))| 



= sup \{iu'\t)-V)Pj,<l),v'\t))\ 

Il0llv=l 

< ||n^(t)||(^.o(Tr))3||P^0||y||i;^(t)||^,. 
Using again fHOj) and the embedding theorem if^(T) C L°°{T), we get that 

dv^ 

^^P ll~irWIU°°(o,r;V") < oo. (41) 
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The estimate (jUJ means that {^^^} is uniformly Lipschitz in V. On the 
other side using the estimate f l5^ . is inside a bounded ball of H. 

Hence, the set {v^{t), Vt} is a compact subset of V. Using the Ascoli- 
Arzela theorem, we can extract a subsequence called again v^{t) such that 

v^-^v in C{[0,T];V') 

and V G C{[0,T]]V' . Moreover, using the estimate (HTl) . the limit v is 
Ltpi[0,T];V'). 

Using dUD and ([381)2 we get that 

ll~7r(^)iU°°(o,2^;V) < oo; (42) 

this result and fBD]) allow to use the compactness theorem (Aubin-Lions). 
Therefore we can extract a subsequence, again called such that 

u = (l - aA)-^v in LP{[0,T]; H^-'), 

for some arbitrary e > and p finite. We shall take e < | in order to use 
that if^~'^(T) C L°°{T). Now, we have all the ingredients to pass to the limit 
in the system (l38l) that we are going to write in the weak form: 



{v^{t,x)-v^{s,x),(f){x))dx= / {{u^{r,x)-V)PN(j){x),v^{r))dx dr. 
T Js Jt 

It is easy to pass to the limit on the l.h.s. of the above equality. Let us focus of 
the r.h.s. of above equality: the non linear term ((m^ ■ V)P/v(/>, v^) converges 
weakly in L^{0,T] V), since converges strongly in L^(0,T; [L°°(T)]'^) and 
converges weakly in L^(0,T; if) (due to fl39|) . considering possibly a new 
subsequence). □ 
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